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We propose a new method to probe for variations in the fine structure constant a using clusters 
of galaxies, opening up a window on a new redshift range for such constraints. Hot clusters shine 
in the X-ray mainly due to bremsstrahlung, while they leave an imprint on the CMB frequency 
spectrum through the Sunyaev-Zel'dovich effect. These two physical processes can be characterized 
by the integrated Comptonization parameter YszD\ and its X-ray counterpart, the Yx parameter. 
The ratio of these two quantities is expected to be constant from numerical simulations and current 
observations. We show that this fact can be exploited to constrain q, as the ratio of the two 
parameters depends on the fine structure constant as oc o?'^ . We determine current constraints 
from a combination of Planck SZ and XMM-Newton data, testing different models of variation 
of Q. When fitting for a constant value of a, we find that current constraints are at the 1% 
level, comparable with current CMB constraints. We discuss strategies for further improving these 
constraints by almost an order of magnitude. 

PACS numbers: 98.80.-k 
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Introduction 

Constraints on the fine structure constant arc currently 
derived from a number of different observations, rang- 
ing from laboratory to astrophysical measurements (see 
e.g.[l|). E.g. CMB data from the WMAP7 satellite in 
combination with ACT and SPT data constrain a at 
- 1% level (68% c.l.) at 2 - 1000 (see also [1, fl). 
Other probes, such as the measurement of the CMB tem- 
perature through the Sunyaev-Zel'dovich (SZ) effect 
in clusters of galaxies Q (in a model dependent way) or 
21cm absorption Q can potentially probe variations of a 
at lower redshifts than CMB. Opening up new redshift 
ranges is useful as theory is not a reliable guide to the 
expected nature of variations in fundamental constants, 
so that variations might be e.g. non-monotonic with z. 
In this paper, we propose to constrain the fine structure 
constant combining measurements of the SZ effect with 
the measurement of X-ray emission in clusters of galaxies. 

CMB experiments such as Planck Q, SPT Q] or ACT 
are in fact currently detecting many hundreds of clus- 
ters through the SZ effect. Some of these are known clus- 
ters, while others are newly discovered, and have been or 
will soon be observed in follow-up campaigns by other 
observatories, such as the Chandra or the XMM-Newton 
telescopes in the X-rays. Thus, measurements of both 
the SZ effect and of the X-ray emission of hundreds of 
clusters will soon be available up to a redshift of 2; ^ 1 

The SZ effect is often expressed in terms of the in- 
tegrated Compton parameter YszD\ (see Sect. II Al for 
a detailed definition), while the X-ray emission of hot 
clusters [kT > 2KeV), mainly due to bremsstrahlung, 
can be characterized by the parameter Yx = MgTx (see 
Sect. II B[) . Both YszD\ and Yx are approximations of 
the thermal energy contained in the clusters, and are thus 
expected to strongly correlate with total mass, weakly 
depending on its dynamical state [l^ - [T^ . In the limit 



where gravity completely dominates cluster formation, 
YszD\ and Yx are expected to scale in the same way 
with mass and redshift as power-laws. Indeed, numeri- 
cal simulations suggest they both have equivalent scal- 
ing relations that are close to be self-similar. Thus, the 
YszD\ — Yx relation is expected to be, on average, con- 



stant at all z |13l. ll5Hl9l |. Furthermore, the same sim- 
ulations show that the scatter on the relation between 
YszD\ and Yx is small, at < 15% level. 

So far, the data are consistent with these predictions. 
In particular, the Planck collaboration found no devia- 
tion from a self-similar behaviour of the YszD\ — Yx 
relation using Planck SZ data and XMM-Netwon X-ray 
data for 62 clusters in the redshift range < z < 0.4 
[20| . Compatible results are also found by [2l|, analyz- 
ing a subsample of 28 clusters from the Planck SZ catalog 
observed in the X-ray by the Chandra telescope. Addi- 
tionally, the SPT collaboration confirmed these results 
using SZ data from the SPT telescope and Chandra X- 
ray data for 14 clusters Q. These analyses do not find 
evolution of the relation with redshift. ^ 

In this paper, we propose to use the observed linear 
relation between Yx and YszD\ to constrain the fine 
structure constant. In fact, the Ysz and the Yx param- 
eters have different dependencies on the fine structure 
constant, so that their ratio strongly depends on a. The 
fact that no deviation from a constant have yet been 
observed in the YszD\ — Yx relation can be used to 
constrain variations in a. 

Currently, the most powerful probe of a in the redshift 
range 0.3 ^2^4 are atomic absorption lines in quasar 
spectra. Tantalizing hints of variation of a have actu- 



^ Rozo et al. [3^ do find a hint of evolution with 2, at the So- 
level, in the scaling relations found by the Planck collaboration 
but attribute this effect not to physics but to systematic effects 
in the analysis of the X-ray data. 
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ally been found in this kind of data. A first analysis of 
143 absorption systems over redshift range 0.3 < z < 4.2 
observed by the Keck telescope suggested a time varia- 
tion of the fine structure constant at the 5-sigma level 
(23 |. a result that was however questioned by an inde- 
pendent analysis (see e.g. [l^). A more recent result 
[25| , that combines the Keck data with a new sample 
of 154 measurements from the VLT suggests on the other 
hand a spatial variation in a. This spatial variation fits 
a dipole with declination (—58 ±9)° and right ascenscion 
(17.5 ± 0.9)h. The significance of this claim is 4.2tT, and 
no systematics that could justify such a variation have 
yet been found. This clearly stimulates the attempts to 
find additional hints of variation of a, preferably with 
probes alternative to quasars. Clusters can explore red- 
shifts lower than quasars, (z ^ 1 versus z < 4), but 
these still correspond to a large range of look-back time 
{ct{z = 1) = 7.7Gly versus {ct{z = 1) = U.lGly)). Fur- 
thermore, current quasar data provide tight contraints 
on the variation of the fine structure constant at redshifts 
higher than > 0.3 [l^l, while clusters can already provide 
constraints at lower redshifts. Potentially, clusters could 
be a better probe at low z, as there are many thousands 
of clusters at these redshifts, where it is harder to find 
quasar absorption systems. It is thus worth investigating 
how much clusters can contribute in constraining a. 

Wc describe the dependence of the SZ effect on the fine 
structure constant in Sec. II Al while the dependence of 
the X-ray data is described in Sec. II Bl The combination 
of the X-ray and SZ effect is described in ll CI In SecHIlwe 
derive constraints from a set of clusters observed, for the 
SZ effect, by the Planck satellite and in the X-rays by the 
XMM-Newton telescope. We present constraints both in 
case a is assumed to be the same for the whole sample 
of clusters, and in case a is assumed to vary with space 
position and look-back time. We conclude in Sec. IIIII 



I. METHOD 

A. Sunyaev-Zel'dovich effect 

Over 80% of the baryonic content of clusters of galax- 
ies is expected to be under the form of intergalactic hot 
gas at temperatures of order T 10'' - 10^ K [l^. The 
ionized gas can inverse Compton scatter CMB photons, 
leaving a signature, the Sunyaev-Zel'dovich effect, in the 
CMB spectrum [l3|- This spectral distortion is propor- 
tional to the Compton parameter y, that quantifies the 
gas pressure of the intraclustcr medium integrated along 
the line of sight: 



n,{r)T{r)dl 



ne{r)T{r) 
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where rrie is the electron mass, Ue is the electron number 
density at distance r from the center of the cluster, T is 
the temperature of the gas, b is the projected distance 
from the center of the cluster, i?^ is the radial extent of 
the cluster and ax is the Thompson cross section, which 
depends on the fine structure constant as 
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(3) 



Integrating the Compton parameter over the angular 
extent of the cluster provides the integrated Compton 
cylindrical parameter y, 



cyl 
SZ ■ 



Y^y\R)D\ = J^' y{b)2Mh ^ D\ y(0)2^0d0 

(4) 

where we expressed the projected distance as 6 = ODa, 
with 6 is the angular dimension of the cluster and Da 
the angular diametre distance, defined as 



Da{z) = (l + z)- 



H{z) 
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Here, a;„i = ^mh"^ is the physical matter density and 
CJA = ^IaIt-^ is the physical dark energy density. In the 
following , we will assume a flat ACDM cosmology. 

The cylindrical parameter can then be linked (see e.g. 
[29} ) to the commonly used spherical integrated Compton 
parameter, defined as: 



rrieC^ Jo 



n{r)T{r)4:Trr'^dr (7) 
(8) 



We will use the spherical integrated Compton parameter 
in the following, calling it simply Ysz- 

Thus the Ysz parameter depends on the fine structure 
constant via the Thompson cross section in Eq. [3] as 



Ysz oc a 



B. X-rays 



(9) 



At the high temperatures of galaxy clusters, 
the intergalactic gas emits mainly through thermal 
bremsstrahlung, whose emissivity, i.e. the emitt ed p ower 
per unit volume at frequency is (see e.g. |26l - [28| ): 
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where Ui is the ion density, k the Bolzmann con- 
stant, h the Planck constant, Z is the atomic number. 
gff{Z, T, v) is the Gaunt factor which corrects for quan- 
tum mechanical effects and varies slowly with frequency 
and temperature. 

We do not observe the emissivity directly, but the sur- 
face brightness, i.e. its integral along the line of sight at 
different angular distances from the center of the cluster 



Yx relation and a 



e^{r)2rdr 



(11) 



Thus, deprojecting the surface brightness allows mea- 
suring the emissivity, which depends on the estimate of 
the angular diameter distance as e^, oc D^^. 

The magnitude of the X-ray emission is often quanti- 
fied through the Yx parameter 13, 29], which is analo- 
gous to the SZ parameter Ysz- It is defined as 



Yx = Mg{R)Tx 



(12) 



Here, Mg{R) is the X-ray determined gas mass within 
a certain cluster radius R, and Tx is the spectroscopically 
determined X-ray temperature of the cluster, determined 
within a cylindrical annulus. The gas mass is defined as 



(13) 
(14) 



Mg{R) = fiemp I UedV 
oc n„R^ 



with fi is the mean molecular weight of eletrons. The 
gas mass can be determined from X-ray data as the den- 
sity profile of the cluster can be inferred from the emissiv- 
ity in Eq. 1101 which is in turn obtained from the observed 
surface brightness in Eq. 1111 This assumes that the tem- 
perature can be spectroscopically determined and that 
the angular diameter distance is known. 

We can then link the inferred gas mass to the fine struc- 
ture constant from Eq. [10] and [14] 



Mg{R) cx Ve^Ti/2e'«>'/fcT 0,-3^3 (^5) 

(16) 
(17) 



From Eq. [T2l [8] and [TU the ratio between Ysz and Yx 
depends on the structure of a cluster as 
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The Ysz parameter depends on the gas mass weighted 
temperature, while Yx depends on the X-ray temper- 
ature. Both are approximations of the same physical 
quantity, i.e. the thermal energy of the cluster. It is 
then clear that if clusters were isothermal, the ratio be- 
tween the two would exactely be equal to a constant. 
However, the ratio can still expected to be constant if 
the temperature profile of the clusters is universal. This 
condition is fuUfilled if the evolution of clusters is com- 
pletely dominated by gravity, weakly depending on gas 
physics. In this case, both Ysz and Yx are expected to 
strongly correlate with the mass of the cluster via the 
virial theorem, both with the same dependence on mass 
and redshift [lllil ^ 

Numerical simulations have shown that indeed the two 
parameters have scaling relations with the total mass of 
the cluster that are very close to be self-similar, i.e. that 
Yx,Ysz oc M5/3£:(z)2/3^ with E{z) = {H{z)/Ho). Fur- 
thermore, they have shown that the scaling relation be- 
tween Yx an d Ys zD^i has very small scatter, at the level 
of ~ 15% p^l - ll7| . The relation between the two is also 
expected to be independent of redshift, as their scaling 
relation with mass have the same dependence on cosmol- 
ogy. Finally, the relation seems not to crucially depend 
on the dynamical state of the clusters [29| . 

Based on all these considerations, the ratio between 
Ysz and Yx is expected to be constant 



YszD\ 



Y- 



const 



X 



for clusters at different redshifts or space positions. 

This fact can be exploited to constrain the variation of 
the fine structure constant at different time/space posi- 
tions i, as 



'YszD 
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(21) 



The dependence on the angular diameter distance de- 
rives from the dependence of i? = 6Da and from the fact 
that the emissivity is derived from the observed surface 
brightness, oc . 

The Yx parameter thus depends on the fine structure 
constant as 

Yx oc a~^-^ (18) 



^ This can be easily shown as follows. Assuming spherical sim- 
metry, the mass M/\^ included in a radius r/^^ within which 
the mean density is Ac times the critical density at redshift z, 
pc(z)= pcmH(z)/HoY, is: 

For the virial theorem, kT oc oc M^^/r^^^, with ip the grav- 
itational potential of the cluster. Thus Ygz>Yx ~ Af^^T oc 
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where 



(YszDl\ 



is a reference value of the ratio that as- 



sumes a reference value of the fine structure constant ag. 
The method enables us to measure the relative variation 
of a with respect to ao as a function of redshift and space 
position. Alternatevely, if one could reliably estimate a 

reference value of ^ ^^Yx^ ) knowing the value of ao, 
e.g. from simulations, it would also be possible to have 
an absolute measure of a for each cluster. 

In any case, if a variation is detected, it could be clearly 
either due to an uncorrected astrophysical or instrumen- 
tal systematic error or due to an actual change in a. But 
if no variation is detected a limit on the variation of a can 
be extracted. We cannot logically exclude the possibility 

that intrisic changes of the ratio ot uncorrected 

systematics might provoke a variation in the Ysz — Yx 
relation that conspires to precisely cancel a true variation 
in a resulting in no apparent variation. This case would 
lead to a false rejection of the variation hypothesis. 



II. CONSTRAINTS FROM CURRENT DATA 
A. Data 

We present in this section constraints on a from cur- 
rent data. For the analysis, we use SZ and X-ray data 
from a subsample of the Planck Earl y Su nyaev-Zel'dovich 
cluster sample [s^, as reported in [2Cl| . The clusters of 
the ESZ sample are detected in the Planck all-sky maps 
through their thermal SZ imprint on the CMB. They are 
characterized by a S/N higher than 6, and are required 
to have a X-ray counterpart in the MCXC catalog jsij . 
The subsample then reported in [2^ is composed by 62 
clusters that had been observed by the XMM-Newton 
telescope, that are not contamined by flares and whose 
morphology is regular enough that spherical symmetry 
can be assumed. We additionally exclude from the anal- 
ysis cluster A2034, whose redshift estimate is discordant 
in and [IH, as noted by [13]. We thus use 61 clus- 
ters, in the redshift range 0.044 < z < 0.44. The Ysz 
and Yx parameters we use here are measured within a 
radius i?500j i-e. the radius at which the mean matter 
density of the cluster is 500 times larger than the critical 
density at the redshift of the cluster. Furthermore, X-ray 
temperatures are defined within a cylindrical annulus of 
radius [0.15 - 0.75]i?5oo- 

Fig. [T] shows the space and redshift distribution of the 
clusters used. 

These data are neither a complete nor a representative 
sample of clusters, and the observation of a larger sam- 
ple of clusters in the X-ray will be required to properly 
characterize the Planck clusters, in particular to study 
the intrisic scatter and Malmquist bias, as well as pos- 
sible systematics [3j]. However, we use this dataset to 



B. Analysis and constraints: Constant a 

We first analyze the data in order to find constraints 
on a under the simple assumption that no evolution in 
time or space is present, i.e. that the YszD\/Yx ratio 
is a constant. Any deviation exceeding statistical error 
is attributed to intrinsic scatter. 

We calculate the mean of YszD\/Yx through a modi- 
fied weighted least square method (MWLS). This method 
differs from a simple weighted least square because it 
takes into account the fact that statistical uncertain- 
ties on YszD\/Yx , calculated by propagating the sta- 
tistical errors on YszD\ and Yx, can be underesti- 
mated or can neglect intrinsic scatter. A weighted least 
square method provides in fact a simple weighted mean 
of \og{YszDA/CxszYx) = -0.056 ± 0.01, with a 
per degrees of freedom of /dof = 233/60. Clearly, 
such a high reduced x'^ might indicate either that a con- 
stant is a poor description of the YszD\/Yx data, or 
the presence of e.g. additional intrinsic scatter. Under 
this second assumption, wc can account for a possible 
wider dispersion of the data by quadratically adding to 
the statistical error of each data point a constant term 
fintr(see e.g. [35| ) for the unknown intrinsic scatter. The 
weighted mean and the intrinsic scatter are then jointly 
determined so that the reduced x^ equals 1. Follow- 
ing this method, we obtain \og{YszDA/CxszYx)i = 
-0.056 ± 0.03, i.c.,{YszDA|CxszYx)^ = 0.96 ± 0.027. 
This result is in perfect agreement with the results found 
by [2^ and [s^- We find that the intrinsic scatter term 
for each data point is equal to (Tintr = 0.18. We note 
here that we do not correct the data for Malmquist bias, 
which for this set of data is not expected to modify the 
best fit [l^l , but might provide a slightly higher estimate 
of the intrinsic scatter compared to corrected data. 

In order to check the results from the MWLS method, 
we also calculate the mean as a simple arithmetic average 
and estimate its uncertainty by bootstrap resampling. In 
this case we obtain {YszDa/CxszYx) = 0.96 ± 0.021, 
in agreement with what previously found. The scatter 
in this case is calculated following [g^: we calculate the 
scattering term a[^^^^ as the quadratic difference between 
the raw scatter draw and the statistical uncertainty 
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(23) 
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where Xr is the reduced , dof is the number of de- 
grees of freedom, in this case equal to dof = 60, and TV 
is the number of clusters, in this case equal to = 61. 



provide a first estimate of the constraints on a that one 
can derive from this dataset. 



The recovered scatter is cr'^ 



0.17 ± 0.026, in per- 



fect agreement with what found with the first method. 
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FIG. 1: Left: Right ascenscion and declination of the Planck ESZ cluster subsample used for the analysis. Right: 
YszD\/CxszYx of the clusters in function of redshift. The error bars shown are calculated from error propagation of the 
errors on Ysz and Yx as published in [20l ]. 



The uncertainty on the scatter is calculated as in [36|, 
A{aU) = al,,{2N[N - 1))-' E(l + T{x.)/'Jintr^? ■ 

We can then convert these results to a measurement on 
a. The current assumption is that {YszD\/CxszYx) 
ratio is constant, and thus that the fine structure con- 
stant has the same value for all the clusters considered, 
a = aQ. The uncertainty on a is then: 



aia) 1 a{YszDyCxszY: 



X) 



ao 3.5 (YszDl/CxszY: 



0.0077, 



i.e. a constraint on a at 0.8% level at 68% c.l. This 
uncertainty includes both statistical error and intrinsic 
scatter, but does not include uncertainties on the cosmo- 
logical parameters used to determine the angular diam- 
eter distance, here chosen for our reference cosmology, 
i. e. flat ACDM with Hq = 70Km/s/Mpc, = 0.3, 
and f^A = 0.7. 

The cosmological parameters are not perfectly known, 
and degeneracies with a could limit the constraining 
power of clusters. We now analyze the impact of these 
uncertainties. 

The dependence of the (YszDj^/CxszYx) ratio on the 
angular diameter distance is 



YszDj 
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-0.5 

(26) 



where ref indicates the angular diameter distance cal- 
culated with the reference cosmology, and true indicates 
the unknown true cosmology. 

First, a wrong estimate of the angular diameter dis- 
tance could generate a "fake" evolution with redshift of 
the {YszD\/CxszYx) ratio. Second, the uncertainties 
on the knowledge of Da can affect the errors on a. Still, 
the dependence is weak and current constraints on the 
angular diameter distance are at the level of a few per- 
cent [33, HI]. We thus expect that, at least for current 



data, the uncertainty on cosmological parameters should 
not affect constraints on a. 

In order to quantify this statement, we redetermine 
the uncertainty on a through the Fisher Matrix method- 
ology, marginalizing over cosmological parameters (see 
e.g. [33 )• Together with the fine structure constant, we 
consider as parameters the physical matter density Q,mh^ 
and the Hubble constant Hq. We impose priors on Q,„ih^ 
and on the Hubble constant Hq from a combination of 
CMB data (WMAP7, ACBAR, SPT, ACT), the galaxy 
power spectrum extracted from the SDSS-DR7 luminous 
red galaxy sample and the constraints on Hq from the 
Hubble Space Telescope observations, obtained when a 
ACDM + a model is fitted to the data [2|. The priors 

2.1Km/s/Mpc. 



are crP{n,nh^) = 0.021 and cr^(i?o) 
The Fisher matrix for Y = YszD\/Yx is 
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E 

k=l 



dYk 



dYk 



(27) 



where 9 are the parameters and a{Yk) is the error on 
the k data point, that includes the statistical uncertainty 
and the intrinsic scatter on Y^ previously determined. 
To calculate the derivatives, we assume a fiducial true 
cosmology equal to the reference one, with a ~ aref = 

Adding the priors, the total Fisher Matrix is 



P,,, with P,, 



1 



2% 



(28) 



where af.^ is the prior on parameter ii. 

The resulting error on a from the Fisher Matrix is 
given by {F~^^^)aa, so that a{a)/ao = 0.0086, while we 
recover the previous result when assuming the cosmolog- 
ical parameters perfectly known. Thus we conclude that 
for current data, the uncertainties on the cosmological 
parameters only marginally affect the constraint on a, at 
the level of - 10%. 
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FIG. 2: YszD\/CxszYx in function of rcos(6l) for the 61 
clusters used in the analysis. 



As more and more clusters are found, this might be- 
come a limiting factor for constraints on a from clusters. 
We estimate that cosmological parameter uncertainties 
become the dominant source of uncertainty on a once 
~ 6000 clusters are observed. In this case the constraint 
on a is a{a)/aQ = 0.003. On the other hand, upcom- 
ing data from on-going experiments such as Planck are 
expected to improve the constraints on cosmological pa- 
rameters. In particular (40| forecasts the constraints on 
cosmological parameters from a combination of future 
CMB and weak lensing data for a Planck satellite-like 
experiment and a Euclid satellite-like experiment, for a 
model where a is also allowed to vary. In this case the 
forecast uncertainties for the Hubble parameter and the 
matter density are a^iHo) = 0.34 and crP(w™) = 0.0007. 
With 6000 clusters the constraint on a would improve to 
a{a)lao = 0.001. 



C. Testing the dipole 

We now present results for a fit of the data that allows 
for a variation of a in space and with lookback-time. In 
particular, we consider the same model of variation used 
by [13] to fit the quasar data. This is useful in order to 
compare the constraining power of clusters to quasars, 
as the two sets of data have different redshift and space 
distributions. The model of variation we adopt for a is 

— = Arcos(6l) + l (29) 
where r is the lookback time, 

J H{z'Y 

9 is the angular distance from the best fit dipole and A 
is the amplitude of the effect. Thus, the relation we fit is 



of the type 

(F)i/3-5 ^ {a + hrcosO) 

b = Y^^^-^A (30) 

We fix here the position of the dipole to the best fit found 
by so that we only fit for a and b with the MWLS 
method, again allowing for intrinsic scatter. The results 
on the derived parameters of interest are Yq = 0.96 ± 
0.028 and A = (-5.5±7.9) x 10-^GyL~\ The constraint 
on A is much weaker than the one obtained from quasars, 
namely A = (1.1 ± 0.25) x lO-^GyL'^. 

As in the previous section, we verify that this con- 
straint is not currently affected by uncertainties on cos- 
mological parameters. With the same Fisher Matrix pro- 
cedure previously described, marginalizing over the un- 
certainties on and Hq, we obtain (j{Yo) = 0.031 and 
a{A) = 8.1 X 10~^. Thus, also in this case uncertainties 
on cosmological parameters impact the constraints at the 
~ 10% level. These findings indicate that current clus- 
ter data cannot directly probe the quasar claim, at least 
under the assumption that the variation is an increas- 
ing function of redshift as in Eq. [SO] This is however 
a phenomenological assumption that is not justified by 
any theoretical model. Thus, variations of a could have a 
more complicated evolution with time. Clusters of galax- 
ies could potentially unveil variations at low redshifts. 

III. CONCLUSIONS 

We propose a new method to constrain the fine struc- 
ture constant by using SZ and X-ray observations, open- 
ing a complementary redshift window on a. With 61 
clusters in our data set, no evolution has been detected 
in the scaling relation between the integrated Compton 
parameter YszD\ and the X-ray analogous parameter 
Yx so far. We can take advantage of this fact to put a 
constraint on the fine structure constant in the redshift 
range < z < 0.5 at the 0.8% level. The ratio between 
the two parameters have in fact a strong dependence on 
the fine structure constant, namely {YszD\/Yx) oc a^-^. 

This constraint is not limited by degeneracies with cos- 
mological parameters at this point. It is conservative 
since we assume that the intrinsic X-SZ relation does not 
evolve. If there were a systematic evolution it would be 
absorbed in our estimate of the intrinsic scatter which 
softens the constraint. 

Further improvements to this constraint arc possible. 
Cosmological parameter degeneracies are subdominant 
until 6000 clusters are observed, at which point the con- 
straint would improve to 0.3%. The intrinsic scatter in 
the Ysz/Yx relation could potentially be reduced with 
improved modeling of the physics in the cluster gas, lead- 
ing to a further improvement. This is currently the 
dominant factor limiting the constraints. Also clusters 
are unlikely to improve significantly the constraints on 
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spatial variation of a. Ultimately a combination of in- 
creased numbers of clusters with improved cosmological 
parameters constraints from upcoming large surveys such 
as Planck and Euclid could sharpen the constraint to 
< 0.1%. This forecast takes cosmological parameter de- 
generacies fully into account. 
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